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Abstract. We present results on two-particle angular correlations in p+p and 
Cu+Cu collisions over a broad range of (r],4>). The PHOBOS detector has a uniquely 
large angular coverage for inclusive charged particles. This allows for the study of 
correlations on both long and short range pseudorapidity scales. A complex two- 
dimensional correlation structure emerges which is interpreted in the context of a 
cluster model. The cluster size and its decay width are extracted from the two-particle 
pseudorapidity correlation function. Relative to p+p collisions, Cu+Cu reactions show 
a non-trivial decrease of cluster size with increasing centrality. These results may 
provide insight into the hadronization stage of the hot and dense medium created in 
heavy ion collisions. 



1. Introduction 



Multiparticle correlation analyses have proven to be a powerful tool in exploring the underlying 
mechanism of particle productions in high energy hadronic collisions. Both short- and long- 
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range correlations have been discovered in the past decades [TJ [2] which have been given 
a simple interpretation via the concept of clustering [21 dj. In a cluster model, clusters of 
hadrons are formed first and emitted independently according to some basic scheme. They 
then decay isotropically in their center of mass frame into the observed hadrons. Two-particle 
angular correlations provide detailed information about the cluster properties, e.g. their 
multiplicity ("size") and extent in phase space ("width"). In heavy ion collisions at RHIC, 
the expected formation of a Quark Gluon Plasma (QGP) could lead to a modification of the 
clusters relative to p+p collisions [5j. This study should provide a useful baseline measurement 
for understanding hadronization stage in A+A collisions. 



2. Analysis method 

Covering pseudorapidity range (rj = — ln(tan(#/2))) —3 < r] < 3 over almost full azimuthal 
angle, the PHOBOS Octagon detector [6] is ideally suited for direct study of the angular 
correlations of the particles emitted from clusters. Following a similar approach as in Ref. [2], 
the inclusive two-particle correlation function in (Arj, A</>) space is defined as follows: 



where F n (Ar], A(p) is the foreground distribution determined by taking two-particle pairs from 
the same event and B n (Ar], Ac/>) is the background distribution constructed by randomly 
selecting single particles from two different events with similar vertex position and centrality. 
The track multiplicity, n, is introduced to compensate for the trivial dilution effects from 
uncorrelated particles. R(Ar], Acj)) is defined such that if a heavy ion collision is simply a 
superposition of individual p+p collisions, and thus has the same local correlations, the same 
correlation function should be observed. 



3. Two-particle correlations in p+p collisions 

Fig. [1] shows the two-particle inclusive correlation function in p+p collisions at yfs = 200 GeV 
as a function of Ar] and Acj). To achieve the correlations between primary particles, a set of 
correction procedures has been applied, based on MC simulations. The complex correlation 
structure suggests that the short range correlation is approximately Gaussian in A-q and 
persists over the full A<j) range, becoming broader toward higher A(j>. If clusters are the 
precursors to the final measured hadrons, a high p? cluster would contribute to a narrow peak 
at the near-side (Ac/> near 0°) region of the correlation function in Fig. [IJ while a lower pr 
cluster will contribute to the broader away-side. 

To study these features quantitatively, the two-dimensional (2D) correlation function is 
projected into a one-dimensional (ID) correlation function R(Arj) shown in Fig. [2j It takes a 
form derived in Ref. [I] in an independent cluster emission model: 



R(Ai]) = a 



r(Ary) • 
B(Arj) 



(2) 



where B(Arj) is the background distribution obtained by event-mixing. Here the multiplicity 
dependence of B(Ar]) is not considered since the multiplicity distribution in p+p collisions is 
relatively narrow with cr(n) / < n > of around 0.25. The parameter a = <K ^ff^^ > contains the 
information about the cluster size K and T{Ai]) is a Gaussian function oc exp{— (Arj) 2 / (4S 2 )) 
characterizing the correlation of particles produced by a single cluster, where y/25 corresponds 
to the decay width of the clusters. The effective cluster multiplicity, or "size" is defined to be 
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Figure 1: Final fully corrected two- 
particle correlation function in (A77, Acp) 
for p+p collisions at y/s = 200 GeV. 
A small region near (A77, A0) of (0,0) 
has been removed due to its large 
experimental uncertainty arising from the 
secondary effects. 
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Figure 2: ID two-particle rapidity 
correlation function in (A77) for p+p 
collisions at yfs = 200 GeV together 
with a fit from a cluster model. The 
fit parameters are: 7^^=2.44+0.08 and 
(5=0.66+0.03 with 90% C.L. errors 



i^efj = — <k> \- 1 =< K > + < £ > • Of course, without any knowledge of the distribution 

of K, it is impossible to directly measure the average cluster size < K >. However, by a x 2 
fit of Eq. [2] to the measured two-particle correlation function, an example of which shown in 
Fig. [21 the effective cluster size K e g can be estimated, as well as the cluster decay width y/26. 
A K c f[ of about 2.5 indicates that on average every charged particle is strongly correlated with 
another 1.5 particles, if it is assumed that aj^=0. 

Fig. [3] shows the collision energy dependence of K c g and 5. PHOBOS data at yfs = 
200 GeV and 410 GeV are consistent with the previous UA5 measurements [I] and show a 
significant increase of the cluster size with energy. The cluster decay width is essentially 
constant with collision energy. The observed cluster size exceeds the expectation from 
resonance decay (about 1.7) p] even at very low collision energy [HE], suggesting additional 
sources of short-range correlations are required to account for the experimental results. 

4. Two-particle correlations in Cu+Cu collisions 

In Cu+Cu collisions, a cos(2A<^>) dependence is observed, and attributed to collective flow. 
However, in this work, the A0 dependence is integrated out in order to just study the 
cluster properties in pseudorapidity. As was explained above, the two-particle pseudorapidity 
correlation function in Cu+Cu is fit by Eq.[2[ The resulting effective cluster size as a function of 
participating nucleons (N part ) is shown in Fig.[3]for Cu+Cu collisions at y/s NN = 200 GeV. The 
dashed line indicates the value found in yfs = 200 GeV p+p collisions, and suggests that the 
cluster properties are similar in the two systems. This seems reasonable if clusters are produced 
in the late stages of the collision evolution and mainly reflect the hadronization process. 
However, it is observed that the cluster size decreases with increasing collision centrality. In 
comparing the data with dynamical models, it is found that AMPT gives the same qualitative 
trend as the data, but is systematically lower in magnitude. By contrast, HIJING remains 



PHOBOS p+p 
ISR p+p 
SPS-UA5 p+p 
PYTHIA p+p 
HIJING p+p 



H 1 — I I I I I 1 1 1 1 — I I I I I 1 1 





• 


PHOBOS p+p 




★ 


ISR p+p 


0.4 


■ 


SPS-UA5 p+p 


O 


PYTHIA p+p 




□ 

i i 


HIJING p+p 
1 



10 2 10 

\S (GeV) 

Figure 3: K e s (top) and 5 (bottom) 
as a function of y/s measured by 
PHOBOS (solid circle), UA5 [lj 
(solid square) and ISR [2l [7] (solid 
star) experiments in p+p and p+p 
collisions. 
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Figure 4: Effective cluster size K e R as a 
function of N part in Cu+Cu collisions at 
^/s NN = 200 GeV (solid circle). The dashed 
line shows the corresponding value in p+p 
collisions measured by PHOBOS. 



constant with increasing centrality. Further comparison of Cu+Cu and Au+Au data with 
p+p should provide more information on the dynamical origins of the centrality dependence 
seen in the Cu+Cu data. 

5. Conclusion 

In conclusion, the two-particle correlation function for inclusive charged particles has been 
studied over broad range in r\ and <p. In particular, it is shown that the observed short- 
range correlations in pseudorapidity have a natural interpretation in terms of clusters. In this 
approach, multiple particles are understood to be emitted close together in phase space, with 
a typical cluster size of 2-3 in p+p collisions. In Cu+Cu, clusters have a similar size but show 
a non-trivial decrease in size with increasing centrality. Future work will extend this study 
by providing a comprehensive study of two-particle correlations in p+p, d+Au, Cu+Cu and 
Au+Au reactions. 
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